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The Primes Lattice 


PROBLEM 63 


The figure on the cover shows a lattice in which the 
heavy lines are the coordinate lines for the odd primes 


foe, 13, 17, 19, 23,29). 


Most of the areas 


enclosed by the heavy lines are rectangles, but some areas 
are squares. In the table given here: 


5 
i 


13 
25 
Ad 
61 
109 
137 
allyl 
253 


25 

a9 
121 
169 
289 
361 
Beg 
841 


520000 
.510204 
2338843 
. 360947 
Neos 
~ 379501 
-410208 
. 300832 


the first column shows the limit of the lattice to be 

considered; the second column gives the area within that 
limit that is enclosed by squares; the third column shows 
the total area up to the limit; and the 
shows the ratio of the area enclosed by 
up to the limit 29, 
squares, 9 4x4 squares, and one each of 3x3 and 6x6, for 
a total of 253 square units out of 841. 


total area. Thus, 


As the limits are increased, will 


fourth column 
squares to the 
there are 16 2x2 


the ratio stabilize? 
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@ Many intriguing problems in numbers involve the 

testing of a given number for primality. The most 
straightforward method for determining whether or not 
a number N is prime is to divide N (odd) by every odd 
prime less than or equal to the square root of N. 


2019-3 


D. H. Lehmer, in a letter to Thomas R. Parkin 
(May 18, 1965) suggested the following improvement: 


"When a number N is being tested for primality and 
fails to be divisible by small primes, say less than L, 
then give up the search which may cost on the order of YN 
units of effort and apply the exponential test 


2N ¢ 2 mod N (A) 


(which implies that N is composite) whose cost is only of 
the order of log N. In case (A) fails to hold, we don't 
know for sure that N is a prime, but then you can always 
pick up the old program at L and continue as before. 


Searching for Primes 


"Come to think of it, perhaps one should replace L 
by I+d; it would be a shame to quit at L and have the 
factor just around the corner!" 

¢ For most primes, 2N = 2 mod N, but there are 
infinitely many composite numbers for which the relation 
holds, although such numbers are relatively rare. For 
example, 


2061 = 2 mod 561 


and other examples are 645, 4371, 10585, 88561, and 137149. 
Oliver Gross has shown that if p is a prime greater than 3, 
then 


LO Ren ame ee 
N= (2 aes 


ys such a number. The generation of the first 7 
mbers is shown here. 


il at ‘ P 
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Lehmer's scheme is one of many ingenious devices to 
make tractable the problem of determining primality of 
numbers that are large (in the sense of lying beyond the & 
range of published tables. Systematic tables of the 
primes up to 104,000,000, and scattered tables in higher 
ranges, are available.) But even so, the value of 
"large" soon becomes excessive. For a number like 


N = 2722258935367507707 7069968594 54145691647 


Asal 
(which is 2 3 1), using the straightforward method of 
dividing by every prime less than or equal to the square 
root of N, and assuming that 


1) all the base primes are known and available 


(there are some 2.6 : 107° of them), and 


2) division of N by any of the base primes can 
be accomplished in one microsecond, 


then the task of verifying a number of the size of N as 
prime or composite could take 85000 years. Clearly, for 
numbers larger than this N, the task of determining 
primality is fairly hopeless without special techniques. 


Just such a technique is furnished by the Lucas-Lehmer 
test for numbers of the form ) 


heat 2) ae 


the so-called Mersenne numbers. In 1644, Father Marin 
Mersenne conjectured that M is prime for only four values 
of p above 19; namely, for p = 31, 61, 127, and 257. For 
p = 257, it was soon established that M is composite. Up 
to 1964, the values of p for which M is known to be prime 
are these: 


2 19 521 4253 

3 31 607 4423 (Selfridge, 1961) 
5 61 1279 eT 

7 89 2203 9941 ? (Gillies, 1963) 
13 107 2281 11213 ) 
17 127 3217 


The values to the right of the line were all verified by 
computer, starting in 1952. (As late as 1950, attempts 


to find a pol me gnTeer than 2+27 - 1 ended in failure. 


The number 3-21 + 1 had been suggested as a likely 

candidate, but it turned out to be composite. The report 

of the research establishing that fact concludes "Thus 

another attempt to discover a larger prime...ends in & 
disappointment." ) B 


The Lucas-Lehmer test proceeds as follows. 


Start 


with A, = 4, Calculate 
2 
A, = [4 a) - 2] mod M 
Mis prime if and only if “Seal =O} The complete tests 
for p = 11 and 13 are shown here: 
M = pit - 1 Me ais - 
Step {Square - 2} Mod M Step {Square - 2} Mod M 
1 4 Il 4 
2 14 14 2 14 14 
3 194 194 3 194 194 
4 37634 788 4 37634 4870 
5 620942 701 5 23716898 3953 
6 491399 119 6 15626207 5970 
T 14159 1877 Tt 35640898 1857 
8 3523127 240 8 Z4U48447 36 
9 57598 282 9 1294 1294 
10 79522 1736 10 1674434 34.70 
git 12040898 128 
12 16382 6) 
not prime prime 


The Lucas-Lehmer test is cheap to perform on a 
computer, compared to any known test for arbitrary 


numbers. Nevertheless, the amount of computation 
increases rapidly with p, and Gillies' results stood 
as the records for 8 years. Bryant Tuckerman, of IBM's 


Thomas J. Watson Research Center, established the number 
M = 219937 _ 1 as the 24th Mersenne prime on March 4, 1971, 


using a 360/91 and 40 minutes of CPU time. 


His description 


of this work is found in the Proceedings of the National 
Academy of Science, Volume 68, No. 10, October 1971. A 
printout of the decimal expansion of the largest known 
prime is shown. Mr. Tuckerman has applied the test to 


all values of p less than 20000. 


The bulk of the work in the Lucas-Lehmer test lies 
“in the squaring of a p-bit number (in the case of the 24th 
Mersenne number, the worst case is the squaring of a number 
of 19937 bits). The reduction modulo M (since M is 
always an unbroken string of 1 bits) can be done by 


subtracting and adding single l's. For example, 
test for M,,, stage 5 produces the number 620942, 
reduction modulo M = 2047 is as follows in binary: 


in the 
The 


% 


- Text continued on page 7 
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THE 24TH MERSENNE 


4315424 


45269 
B4632 
61599 
AALT2 
14958 
23016 
67503 
43417 
$5402 
95410 
AS52646 
59940 
GTO14 
BIG6TA 
71419 
57574 
AR AR? 
06170 
22985 
17947 
65577 
07769 
19944 
QAL74 
A412 
53615 
57781 
?3ATT 
9n191 
12444 
QNtay 
RPA} 
GAAAT 
21291 
ABAKN 
4n911 
SAGAR 
R13Aa0 
OAG2T 
23597 
997H3 
55174 
464301 
49128 
TRO65 
20RK6 
96139 
AA282 
54A15 
47972 
49155 
20331 
03232 
45751 
52950 
RA2A5 
09694 
29930 
07535 


79738 
93825 
76980 
42366 
97324 
00344 
97024 
30896 
96551 
43393 
A4TAS 
87278 
23742 
R7176 
16256 
72625 
45228 
37748 
045746 
02533 
“04n8 
27273 
A4“998 
72330 
99940 
RAGIA 
51681 
AL250 
39818 
29983 
99118 
51955 
61231 
55929 
64174 
43133 
97836 
Baie 
3692 
R4610 
93540 
47501 
05064 
00091 
06433 
66161 
47768 
4B2A5 
00498 
79606 
47663 
51456 
14850 
31567 
39372 
75475 
93480 
87305 
05702 
70621 


81626 
76488 
30006 
75486 
23108 
05435 
01260 
13749 
47057 
39190 
88553 
$2745 
13629 
34948 
08421 
81763 
93476 
61419 
91135 
16515 
89134 
91329 
73548 
B4 238 
91043 
o74a7 
27398 
22279 
21290 
78991 
70210 
10437 
68779 
B8726 
77109 
344Aa7 
53298 
28389 
25338 
5670) 
41470 
35509 
4829 
59169 
7653 
93252 
04260 
28349 
64880 
31695 
12294 
90087 
24875 
91228 
70939 
80582 
06778 
471496 
TW172 
21719 


DOODONONNN 


1 
PRIME IS 2 


48055 23551 
35317 90221 
67281 04984 
$3748 02419 
86265 08132 
36903 89862 
24018 70287 
05766 34090 
54814 59200 
84436 89210 
86633 90254 
10029 63059 
89872 06114 
05549 99616 
28233 80168 
19131 35964 
52454 63153 
62070 41548 
34918 78748 
11431 27880 
98356 985662 
42427 75293 
36042 37173 
01371 10927 
08409 6t741 
06479 56866 
27407 60138 
68993 02674 
B7394 44985 
47125 15365 
$4269 45650 
22546 51574 
47156 13405 
aS aay eae 5 2 
02480 63349 
50679 035H2 
31206 67082 
98040 19973 
22099 09466 
36773 97949 
84870 98546 
49427 60552 
03920 73167 
85856 24217 
99959 94343 
26941 20268 
90308 75482 
70455 71062 
96114 48692 
17658 74142 
67113 48158 
29195 9531159 
63303 11801 
94194 86240 
00970 72378 
05625 95660 
94561 69602 
30121 61656 
86252 84319 
82521 92978 


9) 


9937 


oh 


63379 
73345 
09545 
64350 
16264 
89406 
46650 
52895 
85947 
18424 
62899 
14696 
31410 
30515 
62533 
29547 
57404 
19637 
88923 
20090 
42806 
49082 
46602 
44948 
44708 
90485 
89963 
46255 
91827 
20033 
61384 
09678 
719324 
96641 
10917 
83828 
17149 
90780 
46959 
32024 
54265 
59158 
11307 
65636 
30811 
31159 
38033 
75770 
58603 
01597 
§294) 
94462 
88056 
59403 
27101 
08177 
49439 
86750 
74133 
86978 


00 


19839 
84413 
15881 
29593 
6945) 
17928 
03344 
72290 
26148 
84467 
62132 
38757 
40214 
41225 
45064 
76357 
68786 
93651 
42934 
56808 
32311 
60058 
57694 
35578 
43632 
74782 
86203 
17395 
49999 
60579 
91937 
29084 
95455 
61550 
04827 
35626 
31677 
61463 
19301 
20299 
a2Tt2 
16399 
69513 
24771 
84701 
50894 
45446 
146217 
47371 
9B754 
88164 
94449 
85073 
90744 
24585 
43007 
28137 
73574 
12307 
69167 


MN) 


DONO 


05393 
90952 


1 
5 


6077 
1686 


07770 
72011 
57745 
01603 
36213 


u 


0427 


64328 


1 


$625 


72389 
40346 
31254 
63678 
22894 
23192 
93401 


4 
9 


5650 
8520 


58847 
34723 


1 


$037 


42794 
65521 


L 
y, 


8776 
5424 


00722 


9 
el 
3 


35908 
GI43 
3602 


09528 


1 


4060 


64122 


8 
3 
6 
6 
tS) 


6083 
5643 
7541 
1626 
7077 


05730 
80523 


1 


8922 


32898 


4 
o) 


7158 
6715 


56630 
61565 
43636 
27342 
32825 
32352 
56984 
$2321 
38376 


1 


9174 


63495 
95558 
61788 
34625 


NN) 


WS 


doognoanonn 


— 


50432 
82691 
13296 
62363 
78122 
12083 
70315 
60005 
87555 
27684 
24257 
73553 
98090 
52970 
03467 
37019 
58856 
81799 
45170 
68188 
43682 
66523 
55063 
58684 
76892 
10547 
303646 
69831 
35924 
60167 
24503 
59382 
o525t 
91607 
82811 
23657 
40379 
56710 
09792 
31017 
09003 
19067 
96593 
16785 
71281 
13845 
56192 
26272 
59194 
80267 
044646 
48817 
15805 
67801 
78338 
68126 
65759 
82340 
67975 
61843 


26711 
54609 
98437 
90904 
82829 
36148 
43129 
71630 
71168 
60185 
48035 
44759 
96281 
07721 
08061 
34983 
54608 
90066 
57171 
717266 
63294 
15095 
01744 
96465 
05620 
40113 
87627 
63786 
23336 
Ba07G 
96626 
25780 
80101 
98322 
86632 
54820 
28972 
18463 
osca? 
69258 
33847 
96107 
50901 
48246 
61493 
19588 
41308 
51534 
01396 
23452 
66927 
59990 
18118 
93818 
16102 
55955 
98474 
30398 
06784 
01754 


50516 
14801 
42134 
T94GR4 
444TT 
OR44T 
29960 
87519 
44445 
29149 
TRH24 
T9734 
AQ15B8 
146273 
Apa 
517A 
54205 
217THA 
4AlR11 
6N9R4 
15290 
74170 
11LAAT 
AST49 
04272 
94592 
Antsy 
30001 
SATRS 
97546 
77999 
TANRR 
220A7 
22298 
44590 
AR4TI 
43AGI3AD 
406467 
47175 
2NGh? 
9n833 
R4993 
86230 
BATA 
9442) 
AIT 
A7445 
07407 
27063 
4Al24 
Al1a7 
71191 
71N77 
W219 
43975 
M2174 
R5773 
QRT4LA 
246019 
54903 


bub 


$2505 
90078 
62179 
47692 
50226 
48291 
25187 
13739 
TA9TS 
25270 
35806 
46310 
90645 
13557 
54327 
62164 
ANG24 
46716 
25795 
16348 
75297 
77831 
414612 
91777 
27961 
96221 
SR464 
iene) 
N6711 
pazai 
38402 
02738 
25877 
27240 
71457 
AI9DS4L 
R6T44 
37026 
67020 
10702 
42506 
58089 
94810 
62495 
2Th61 
47982 
27546 
62540 
66851 
35691 
47449 
35755 
86539 
09215 
86856 
76709 
$5399 
Fess 
54367 
86411 


41403 
75343 
03963 
31397 
80488 
35473 
78079 
79555 
08862 
28092 
08154 
06783 
69393 
a4o0al 
07144 
29496 
68987 
71134 
88888 
83684 
27983 
91001 
92438 
28699 
63866 
77502 
00427 
22715 
86568 
78570 
38678 
55261 
8974) 
614678 
87138 
63838 
13989 
03777 
65051 
21251 
04119 
66832 
SISSY | 
30251 
42282 
74887 
68143 
51499 
3by29 
57307 
48983 
93338 
‘Soh 
54607 
89423 
22460 
09575 
45506 
60305 
15854 


33068 
74139 
9134) 
83013 
05782 
28367 
01193 
04746 
77996 
26975 
69654 
67393 
44833 
49309 
78876 
O7571 
37243 
71632 
92774 
90562 
434629 
61684 
95814 
26744 
91498 
57556 
98580 
14060 
839'8 
952555) 
16864 
Splat 
15817 
32528 
35123 
195 32. 
18554 
65347 
39542 
41204 
50303 
99297 
95603 
36036 
62909 
92618 
01548 
319AG 
96928 
AMEN) 
B5064 
21217 
29601 
68444 
06960 
86674 
Siete} 
09577 
99034 
29504 


01376 
62968 
28520 
77820 
00287 
27787 
75902 
81543 
48730 
38426 
69325 
32740 
30994 
86630 
86186 
90918 
69214 
71548 
95426 
18222 
44650 
75685 
15491 
21836 
05489 
58110 
69186 
34129 
67477 
54130 
09962 
20440 
04826 
92479 
6B422 
17452 
16612 
A1733 
36075 
29322 
ooo17 
68126 
05240 
34127 
95005 
51417 
77108 
49445 
6949) 
37397 
37578 
06191 
43729 
SUE BIL 
91540 
77445 
20080 
36921 
07084 
56992 


58091 
01920 
56276 
78571 
64659 
95656 
86317 
33253 
43084 
77025 
63833 
21499 
16963 
65733 
19833 
05462 
45092 
17958 
97714 
62933 
99922 
65867 
00609 
96211 
83112 
67452 
26402 
90218 
04960 
49029 
01587 
75620 
55889 
71651 
61380 
25026 
Zoe \sytie) 
67084 
agoel 
53043 
04002 
24423 
7B716 
68366 
57469 
81997 
7772A 
91064 
80517 
81041 
75073 
47718 
40865 
85595 
24998 
20875 
90408 
55919 
81464 
09056 


Bryant Tuckerman 


Thomas J. Watson Research Cer 


March 4, 1971 


¥ 


i 


13045 
11448 
19600 
24190 
39916 
48307 
10841 
47499 
50484 
Leb) 
32767 
30968 
22958 
67719 
20777 
51141 
31537 
77005 
99145 
98654 
06368 
31928 
75221 
Soi 
12446 
20144 
68612 
15706 
01427 
27192 
90905 
32678 
71438 
99369 
07462 
82372 
93566 
A4473 
59835 
17896 
88789 
14008 
38092 
45617 
81053 
90344 
01108 
14660 
25568 
62765 
76496 
50549 
27040 
13613 
79074 
60785 
30036 
54808 
60727 
3674) 


oe 


13629 
64809 
51310 
54474 
47662 
84690 
49642 
10462 
22342 
28954 
07698 
77829 
77995 
11728 
29064 
43666 
16984 
38269 
49623 
86456 
78136 
60882 
6BB22 
101048 
76399 
69819 
68618 
59650 
75406 
22018 
86549 
06244 
kone 
89519 
19140 
44136 
68612 
47020 
35754 
16267 
29414 
65703 
19164 
50770 
2066) 
17285 
60043 
53430 
18508 
Bu SiG 
345146 
36632 
70219 
30424 
53461 
90623 
44649 
15514 
89554 
53903 


LID LAD PDL II 


31858 
02661 
66463 
33284 
65200 
91169 
47337 
48132 
06292 
40120 
99439 
67413 
84899 
53987 
48066 
38418 
07168 
43934 
91639 
69080 
71540 
07017 
30887 
32785 
93195 
91968 
0838 
53260 
25331 
41725 
42350 
48803 
44080 
18780 
04818 
32757 
6582? 
56866 
10321 
04777 
04603 
34218 
43375 
Ata 
00015 
59860 
25892 
53785 
29882 
78598 
86253 
21115 
24383 
22061 
31192 
34750 
22194 
03591 
95487 
09680 


HNONOOEEODORNANNNOONEEARNAANNNONOENDS 


LOOAOtALLOGL LOO oO lato = 620942), 4 


Si +1 
Mor miewato oO Seen. 
2 Se eee Subtract M by 


ae erie a oe 0 1 Oral On subtracting (M+1) 
1 1 and adding 1. 
0) 
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1 


© em 


LOlOL® Aa aa 
iL i 


tee eran ee emer. (ON). 


For small numbers (say, from 9 to 20 digits), there 
are various schemes for determining primality, such as the 
one by Lehmer at the start of this article. For very 
large numbers of a special form, there is the Lucas-Lehmer 
test. In the mid range (20 to 30 digit numbers) there is 
still a need for new techniques. Two of the latest are 
reported in Mathematics of Computation, April 1974: 


D. H. Lehmer and Emma Lehmer, "A New Factorization 
Technique Using Quadratic Forms." 


R. Sherman Lehman, "Factoring Large Integers." 


The latter article reports the factorization 
29742 3156994.06748437 = 372173423 79915202819 
in 122.6 seconds on a CDC 6400, 


PROBLEM 65 


Two Primes Problems 


The 24 odd primes less than 100 are to be placed 
on the 24 faces of four cubes, in such a way that 


1) Any toss of the four dice produces a sum that 
is divisible by 4; or 


2) Any toss of the four dice produces a sum that 
is not divisible by 4. 


Are either of these arrangements possible? ste 
the 24 odd primes are taken to be those from 5 through » 
101, is either arrangement possible? 
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A succession of random numbers (uniformly distributed 
in the range from 001 to 999) is drawn. The numbers are 
progressively totalled until the sum is a prime number, 
at which time the game ends and the score is the number 
of numbers drawn. Problem: what is the distribution 
of the scores? 


PROBLEM 66 


Casual calculation (by hand) indicates that the 
games are short. For 30 games, the distribution was 
as follows: 


Length of game: 


1 BU By Bes Bis reece Shon ae ae) as AUG 
Occurences: pees) eee aL 


1 dl al a a 


Se cca 


Problem 59 (PC18-8) calied for finding 100 digits 
of a number which, when squared, would have its 1.00 
low-order digits the same as the original number, Such 
a number was named automorphic by Maurice Kraitchik, in 
his book Mathematical Recreations. In the 1942 edition 
of that book, the low order digits were given as 


2 
Q 6259918212390625 
= but with an error in the digit indicated. 
iF 
a 
J Sanford Goldfarb indicated a 
straightforward algorithm for 
) extending the number: for any portion 
n) of the number that is known, calculate 
its square, as shown here. The 
Known digits reappear in the product, 
2 and the digit to the left of those 
w that repeat is the next digit sought. 
za 
a : 
Using just this method, a program was written for 
A) an IBM 1620 that produced the following result: 
£ 
a 


3953007 319108 1698029 385098900621665095808638110005 Oo 
574.234.23230896109004 106619977 392256259918212890625 


PC1y-9 


@ 25 Squares 


In the array of 25 squares 
shown at the right, the 
25 primes less than 100 
have been inserted on the 
east side of each square. 
The other 75 numbers from 
OO to 99 are to be placed 
at the north, south, and 
west sides of the squares, 
with two simple 
restrictions: 


(1) The sum of the four numbers in each square is 
not to exceed 210, and 


(2) The sum of the two numbers at each of the 40 
common borders is not to exceed 105. 


PROBLEM 64 


Thus, the problem solution 
could begin with an 
arrangement like this 

one (the numbers have been 
placed facing outward in 
their squares): 
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_ can be extended to chains of 4, 5, 6, or 7, leading to longer and 


The following six prime numbers: 


203425 Chained Primes 
2034239 

ee ne rim 
508559 

254279 PROBLEM 67 
127139 


are chained; that is, for each of the first five, (p-1)/2 is also a 
prime, The chain of 6 is the longest known, and Problem 67 is to find 
a chain of 7 or longer. There is no a priori reason to believe that 
longer chains exist, but it is likely. The simplest such chain is 
eels 58 and 2 


In searching for chained primes, it is not necessary to examine 
all consecutive primes. Table A shows the leading primes for chains 
of 3, which are quite abundant. it can be seen that the differences 
between successive entries are all multiples of 24, and that the 
entries themselves are all of the form 24K+23; it is not difficult to 
show that this is true. Given an entry value, the subsequent 24 
integers can be expressed as the remainders modulo 24; that is, as 
numbers of the form 24K+M, and only those for which M=1, 5, (6, 
13, 17, 19, or 23 could be primes. But since (p-1)/2 and (p-3)/4 
must also be prime, the situation reduces to: 


jo 
(p-1)/2 = 
(p-3) /4 = 


Vay 2 Df Me) 25 © fies 2h gg 
S 6 @&  G il free 12 
2 4u 5 mod 6 


137279 

13999 

145007 

146519 

147047 

148199 

148727 A table of consecutive prime numbers, p, 
249519 

151967 for which (p-1)/2 and (p-3)/4 are 
157679 

159119 also primes. 

165527 

166487 A 
166679 

166919 

167879 

168527 


and the only case that survives the sieve is 24K+23, This reasoning 


longer jumps that can be made among the integers in searching for long 

chains. Thus, for chains of 4, the eligible numbers for the start of 

the chain must be of the form 48K+47; for chains of 7, the eligible 

numbers must be of the form 384K+383. (va) 


ES TD 


More on Problem 43E (PC13-6) and the 44 terms shown in 
PC17-16. The problem, as stated: 


Take the integers from 3 to N. Circle the 3 and 
cross off every third remaining number. Circle 
the next remaining number (4) and cross off every 
third remaining number. Circle the next remaining 
number (5) and cross off every third remaining 
number. List the circled numbers, and give the 
1000th circled .number. 


Sanford Goldfarb pointed out that the 44 terms 
given as a partial solution follow this pattern: 


(3/2)a, - 1/2 for a, odd 


a 
n+i 
even 


aay (372 Jae - 1 Comme 


{but that the 44th term given had a transposition error 
in the last two digits). Thus, the 1000th term would 
be approximately 

Ig 


(3/2) or, antilog 175.9152 


David Ferguson (of Group/3) analyzed the problem 
further, as follows: 


"Suppose instead of taking the integers from 3 to 
N, one were to take the Integers from 4 to N with all 
multiples of 3 omitted. Then clearly one would get the 
same circled numbers in the same order with the integer 3 
missing. Consequently the number of integers eliminated 
preceeding the nth circled number (c,) would be the same 
as the number of integers eliminated preceeding (cy) in 
the original sequence. Hence the number of integers 
eliminated from c,_ 1 to cn in the original sequence is 


exactly the same as the total number of integers preceeding 


Cy, Which were eliminated by 3. But since no c, is 
divisible by 3 (except 3), this result can be stated as: 


G = Cp = [¢)/3] 


solving this for c, there are two cases: 


Problem Solution 
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These last two results can be combined into: 
a 3c es] 
ec, = nel 


@ 


Mr. Ferguson coded that solution for his System/3, yielding 


. the first 1000 terms in 6 seconds of CPU time. The 1000th 


term is exactly 
13344154'704150796868591877 25459711777 1886 35805769080462194.27 
45785807402551367719517 3706881644 37028253 320782974 3715687939 
78924 313240762285 34214982057 14 287 3941052434787 307231 33193 


(177 digits), and the 44th term is 60540698. O 


SSS A AMMOMH 


| 


Among the sieve problems (PC13-6) was that of Ulan, 
which calis for forming the sequence: s 


a 2, 34, 6, Sold, 5135 Loyal oe eee 


in which each new member can be formed in one and only 
one way by adding two different earlier numbers. The 
number 27 will not appear, since it can be formed in two 
ways (26 + 1, 16 + 11); the number 35 will not appear 
since it cannot be formed at all. 


Associate Editor David Babcock found the 1000th 
term of the sequence to be 12294, using an assembly 
langua program of 50 instructions. The first 110 
terms acres the 26) are listed here to use as test data 
for students who wish to try extending the sequence. 


28° 36. 36 4T -48 . 53.) 57. (62). Goenemmean 
82 87 97 99 loz 106 114 126 131 138 145 


Problem Solution 


The Future of Programmers 


What jobs are there that rate the title "programmer"? 
The loose term can include any of the following: 


1. Applications programmers. These are the people 
who write a new tape merge (in COBOL); who computerize the 
company's personnel file procedures; who maintain the 
payroll procedure in its weekly changes; who use Fortran 
to analyze a market survey; and so on. At the lowest 
level, these people are coders and, however skillful they 
may be at coding, that skill is not worth much in the 
marketplace. They include the users of the specialized 
languages (SIMSCRIPT, FORMAC, etc.). 


2. Vendor system programmers. These are the people 
who create new versions of an operating system; who write 
the compilers (i.e., the translators); who create the 
support software for group 1. By definition, they work 
for the makers of the hardware. Their product may be 
bundled or unbundled, but traditionally there seems to be 
no pressure to make it satisfactory to its users. 


3. User system programmers. These are the experts 


who are the interface between group 1 and the vendor-supplied 


software. They rectify the bugs found in the operating 
system and the compilers. They moderate the battles 
between the applications programmers and the vendor's 
representatives over the questions of whether the faulty 
output is hardware or software trouble and, if the latter, 
whose software. They may modify the vendor-supplied 
software to fit their organization's particular needs. 
They may also write systems software when the vendor's 
is inadequate and/or there is none in the marketplace 
that satisfies their organization's real or imagined 
needs. 


4. Software company programmers, These are, or 
should be, the real pros who create the packaged and 
proprietary software for canonical tasks. They differ 
from the group 2 people in that their output must be of 
higher quality; that is, the software they produce must 
be responsive to the user's needs. 


Group 2 programmers will be with us for a long time. 
A few hundred competent ones could probably satisfy the 
needs of the industry for the next decade, but we will 
have several thousand of them. Similarly, group 4 
programmers will be needed for some time but, by 


definition, they are few in number. » 
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» Group 3 programmers will also be needed. Each 
installation will need one or more, particularly during 
the time when vendor-supplied software is as badly written 
and poorly documented as it is in 1974. When support 
software is created that is relatively error-free, human 
engineered, and idiot-proof, this group will dwindle. 


Group 1 is presently the largest group. Applications 
programming is what students (in universities, community 
colleges, and trade schools) have in mind as their eventual 
work. We already have applications programmers by the 
tens of thousands, and we are busy creating thousands more 
each semester. 


The thesis of this article is that we are training 
too many people for blind alley jobs; that the market for 
these programmers is already low and is getting lower. 
Three trends support this thesis: 


1. Packaged programs. In the scientific/engineering 
domain, there are few problems that have not been programmed. 
These programs are far from perfect, and are not available 
in all languages for all machines, but there is little 
demand for a new program to invert a matrix, or to calculate 
standard functions. Even large problem situations (motor 
design, lens design, gear design, etc.) have been packaged, 
and more such situations will be covered. Although it is 
fun (and instructive) to write a new PL/I program for 
solving simultaneous equations by Gaussian elimination, 
there just is no market for the result; we've done that 
task, over and over. 


In the business area, generating programs exist 
for all the stock jobs: sorting, report writing, file 
maintenance. Indeed, through large packages like 
Informatics! Mark IV, nearly every routine task of file 
manipulation can be performed without the necessity of 
writing a single instruction in the usual sense~-and the 
various tasks of file manipulation cover 95% of everything 
business wants to do with computers. 


2. It has been demonstrated (via systems like JOSS) 
that conversational computing can be a powerful tool to 
couple the man who has a problem directly with the machine, 
so that he can be his own programmer, The so-called 
interactive languages (e.g., on-line BASIC, or Fortran) 
are only a first approximation to conversational computing; 
that is, they are largely just remote job entry systems, 
and substitute an expensive terminal for an inexpensive 
keypunch,. To be sure, they cut the turnaround time, but 
they are a long way from being conversational. True » 


» 


conversational computing, which has existed up to now only 
in time-sharing mode, is very expensive, which is why it 
has almost completely disappeared. When conversational 
computing returns, in dedicated systems consisting of a 
mini computer and a terminal, the effective coupling of 
man and machine will again act to eliminate programmers 

as the middlemen. 


3. There are problem situations in abundance for 
which turnkey systems~-~again centered around a mini--will 
appear, The paperwork problems of the operators of small 
businesses (doctors, lawyers, accountants, automobile 
salesmen, real estate salesmen, mailing list users, and 
so on) are sufficiently alike by types to be packaged. 


These three trends add up to a bleak future for 
those who regard programming as a lifetime of work writing 
pretty routines. What the industrial world needs~-pretty 


badly-~-is people who understand computing and data processing, 


and the writing of programs is the smallest part of that 
discipline. Of far greater importance is knowing what 
should be computed; how that computation can be performed 
efficiently and effectively; and--most important--how to 
verify that the computation is correct, If any of those 
elements is disregarded, the end result is not worth much, 
as, for example: 


The correct results, done inefficiently; the 
1 hour run that could have been done in 10 seconds; 


2 The correct results done ineffectively; the 
results are not presented in a meaningful way; 


The correct results done by computer that 
3 could have been done analytically, or by looking 
up previous results; 


Incorrect results, with the programmer unable 
4 to differentiate his computer-produced garbage from 
good work. 


Sadly, the management of our industry seldom asks 
(or can tell) whether work has fallen into one of those 
four categories. ® 
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technological obsolescence among programmers. (A symposium 
in December 1973 on “Exploring the Future" dealt at some 
length with this topic.) The industry tends to force 
programmers into specialization, while raising their salaries 
year after year. Presently, a point is reached at which 
the individual finds himself out of date and replaceable 

by someone half his age at half his salary. That point 

now faces many men at age 50. Partly, the fault lies 

with the system of wages and responsibility we have 
engineered in our industry. Partly, it lies with the 
individual who has not kept up with his trade. In any 
event, it foretells trouble for many people, and again 
points up the dreary future for large numbers of programmers. 


© ; An even more serious problem is the tendency toward 
<a 
a 
Q 
Aa 


If the above analysis is even partly correct, then 
there will have to be serious revisions made in our training 
and education practices. The current emphasis is hundreds 
of courses is on the mechanics of programming (that is, the 
writing of instructions in some higher level language), 
usually on well-defined problems. This training may be 
necessary, but it is far from sufficient to satisfy the 
future needs of the computing world and, at its best, is 
not education. 


---BY FJG 
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19 


Log 19 1.278753600952828961536333475756929317 951129337 39450 
ifgy JUS) 2 9444 389791664404600090274 31887853537 237 37926129913 


VI9 =—- 4. 35889894. 354067 355223698198 385961565913700392523244 
VI9 2. 668401648721944867 339627 3719708 3033509587856918310 
WI9 1.80198312731714230518255395296189025894 370970228005 
VI5 = -:1..52292699821875339338696819152719299862141820457222 
Wi9  1.34237965096210479809378379300029234668987 798576376 
‘VAD —«1.02988216191782398114014848587167216575078284094775 
ig a eee 178482300 .9631872608449100337887 227038836197 33165166 


a es wattas : vii 4261953832015756927 522 39060170044051350752 
an b 


K 2792563049. 59784556528989784012096391704900630586477 
: 32 233098 3964012 328725477 518397 358909651344 


